function [ssd] = SMM()
% calculate the sum of squared difference between data and simulated data
% moments
%   inputs: 
%   outputs: ssd - sum of squared difference of moments

% get g function
g = 

% calculates ssd
ssd = 

end

